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Abstract — Eugene Initial reservoir pressures from well test information were found to be higher than the calculated 
hydrostatic pressures. Three wells in Block 276 PI -sand had higher observed pressures than what was expected, implying 
pore pressures greater than anticipated. It is thought that under compaction in shaly sand and thermal dissipation from salt 
diapirs are the main reasons causing over pressurization. Both effects have been studied and under compaction was found to 
be the dominant mechanism causing pore overpressure. Transgressive shale in PI -sand, resulting from ocean transgression, 
acted as plastic barriers; thereby, negating the transfer of litho static load to deeper formations. This has resulted in constant 
porosities with respect to depth. A developed COMSOL multiphysics model has also indicated that temperature effects from 
underlying salt diapers have also affected pore pressure and that pore over pressurization by conduction was found to range 
from 3-15%. 
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I. Introduction 

EI330 Normal reservoir pressure is the pressure necessary to sustain a column of water to surface (Fertl, 1976]. Normal 
pressures range between 0.43 and 0.50 psi/ft. Normal drilling muds weigh about 9 ppg (pounds per gallon) and exert a 
bottom hole pressure of approximately 0.47 psi/ft of depth. Overpressure refers to pressures higher than normal. Pressures 
lower than normal are called subnormal. Dickey et al. [1968] identified the location of abnormal pressures in southern 
Louisiana and noted that the continental and deltaic facies contain sandy beds and that the neritic (nearshore) facies also 
comprise a few silty and sandy beds that connect laterally to the deltaic facies. The authors added that shelf facies contains 
almost no sandy beds, and the pore fluids cannot escape. It was also highlighted that growth faults are seals that stop the 
lateral flow of pore water toward the neritic facies. 

Abnormal pressure is defined as any departure from normal hydrostatic pressure at a given depth [Bruce and Bowers, 2002]. 
Abnormal subsurface pressures, either overpressure (geopressure) or unde rpressure, are found in hydrocarbon basins 
throughout the world in all lithologies, from all geologic ages, and at all depths [Fertl et al., 1994]. 

Timely and reliable detection of geopressure is imperative to avoid or mitigate potential drilling and safety hazards, e.g.: 

■ Shallow water flow 

■ Blowouts 

■ Shale instability 

During drilling, advanced warning of approaching geopressuring enables mud engineers to adjust the mud weight to prevent 
well and reservoir damage and to determine casing points. This is a particular concern in deepwater wells where the pressure 
difference; i.e., the operating window, between the hydrostatic gradient and the fracture gradient is very narrow. 

For most cases, shale and shaly sands are formations that exhibit over pressurization. Shale is abundant in Block 276 PI -sand 
with an average net to gross of 46 to 48%. A better understanding of porosity mechanics in Block 276 is of interest for many 
reasons, namely: 

1) Shale is the dominant lithology in the block of interest basin. 

2) Expulsion of pore waters from shale could affect the pressure distribution and flow of pore waters and hydrocarbons 
throughout the block of interest basin. 

3) In many cases, shale porosity depends only on effective stress (lithostatic stress minus pore pressure), and excess shale 
porosity provides a useful measure of pore fluid pressure. 

4) Over pressured shales pose drilling hazards and stimulate slumping and fault movement. 

5) Shale porosity profiles can replicate the history of over pressuring and the time of seal formation in the area. 
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To interpret shale porosity, we need to know how shales compact when pore fluid pressure is hydrostatic. Hunt et al. [1998] 
have suggested that, in wells where fluid pressure is hydrostatic, porosity decreases linearly with depth until porosities of 
about 10% are reached. 

At greater depths, as in the case in Block 276 PI -sand, porosity changes slightly. It goes up and down, but does not portray 
signs of degradation with depth. The boundary between linear (stage 1) compaction and no (stage 2) compaction occurs when 
the shale pores have been condensed to the thickness of about 3 mono -layers of water. For shales with a high percentage of 
high surface area clays (illite + smectite + illite- smectite) the limit of normal compaction can be 15-20% or even greater, 
according to Hunt et al. [1998). For shales with a low percentage of high surface area clays, the limit of compaction can be as 
low as 3%. Once shales reach stage 2 compaction, shale porosity is independent of effective stress. If pore pressure is then 
increased (by hydrocarbon maturation, for example), porosity will not increase and overpressures will not be revealed in 
shale porosity. Hunt et al. [1998] showed several examples where the top of overpressure lies in stage 2 compaction zone and 
porosity remains unchanged, despite the large surge in pore pressure. The perception that there is a natural limit to shale 
compaction provides a sound justification of these interpretations. 

This contradicts with a considerable literature suggesting that porosity data can be used to predict excess (above hydrostatic) 
pore pressures [Magara, 1978; Fertl, 1976; Bray and Karig, 1985; Shi and Wang, 1988; Bangs et al., 1990; Luo and Vasseur, 
1992, 1993; Bour et al., 1995; Gordon and Flemings, 1998; and Stump et al., 1998]. 

In addition, Revil et al. [1998] examined porosity-depth profiles and mud weight data in 89 wells in the Eugene Island area to 
determine whether the intervals of approximately constant porosity are caused by compaction reaching its natural limit or by 
fluid overpressure. By comparing the fluid pressures predicted from shale porosity to those computed from mud weight data, 
Revil et al. [1998] found that constant porosity intervals are in a state of disequilibrium compaction. The porosities were 
controlled by pore pressure as Gordon and Flemings [1998] reported, from an analysis of several wells in the area. 

To analyze and infer the pattern of compaction, Revil et al. [1998] developed theoretical relations against which departures 
from hydrostatic compaction can be measured. They found that the pattern of subsurface shale porosity is spatially logical on 
a local scale but varies over the study area. They defined a surface where porosity commences to depart from the hydrostatic 
compaction trend. This surface parallels the H. sellii transgressive shale over the study area. It was proposed that the 
sediments become impermeable when hydrocarbon fluids entered the area and capillary seals formed [Cathles, 2001 and 
Shosa and Cathles, 2001], but alternative explanations are also likely. 

II. Reservoir Simulation and History Matching 

In the attempts to history match reservoir pressure in PI -sand, pressure profiles (Fig. 1) show a significant divergence 
between wells #1, #12 and W-12 observed (solid) and estimated (dashed) pore pressures, during early production time. Early 
departures are indicative of hydrocarbon overpressure caused by under compaction and possibly heat dissipation from salt 
diapers in the vicinity of Block 276. The late time matching could possibly be explained by wells producing under normal 
hydrostatic regime after over pressured hydrocarbons have all been depleted. 



Time (days) 

Figure 1: Pressure profiles in PI sand, observed and simulated 


Page | 42 


International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-3, Issue-7, July- 2017] 


In this work, it is intended to study under compaction in PI -sand, thought to be the dominant mechanism behind overpressure 
in Block 276. Temperature effects on pressure are also investigated. 

III. Arrested Compaction 

In Fig. 2 below, the depth at which porosity departs from normal compaction trend is what is referred to as fluid retention 
depth. At that particular depth, porosity remains constant as is the case in Block 276 PI -sand depths of 7621 to 9600 ft. It is 
also the depth at which pore pressure departs from the hydrostatic pressure trend and fluids become over pressured. The main 
attribute is what is called disequilibrium compaction. 

Porosity (%) Pressure (MPa) 



Figure 2: Relationship between (a) porosity and depth and (b) pressure and depth for a shale that 
becomes over pressured by disequilibrium compaction at the fluid retention depth. Modified from Harrold 

etal. [1999] 

According to Lash and Blood [2006], the fluid retention depth in the Gulf of Mexico varies from as shallow as 0.3 mi. to as 
deep as 0.9 mi (Fig. 3). At that depth, shaly sandstone composition sees a drastic shift in its mineral structure. Smectite 
replaces most of the illite leading into water expulsion. That translates to a higher hydrocarbon storativity and higher pore 
pressure. 



Sedimentation Rate (m/MY) 

Figure 3: Fluid retention depth as a function of sedimentation rate, Lash and Blood [2006] 

Besides, when the sediments subside rapidly as in Fig. 3 above and in Block 276 PI -sand, fluids in the sediments can only be 
partially expelled. The remaining fluid in the pores must sustain all or part of the weight of overly sediments, causing the 
pressure of pore fluids to increase (abnormally high pore pressure). In this case, porosity decreases less rapidly than it should 
be with depth, and formations are said to be under-compacted or in compaction disequilibrium. 
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Furthermore, and according to Hunt [1996], disequilibrium compaction is caused by stepwise transformation to mixed -layer 
clays and illite. This releases more water into the rock and causes a slower compaction of smectite-rich rocks. Overpressures 
can be generated by many mechanism, such as compaction disequilibrium (under-compaction), hydrocarbon generation and 
gas cracking, aquathermal expansion, tectonic compression (lateral stress), mineral transformations (e.g., illitization), and 
osmosis, hydraulic head and hydrocarbon buoyancy [Gutierrez et al., 2006; Swarbrick and Osborne, 1998]. In nearly all 
cases where compaction disequilibrium has been determined to be the primary cause of over pressuring, the age of the rocks 
is geologically young. Examples of areas where compaction disequilibrium is cited as the primary reason of abnormal 
pressure include the U.S. Gulf Coast, Alaska Cook Inlet; Beaufort Sea, Mackenzie Delta, North Sea, Adriatic Sea, Niger 
Delta, Mahakam Delta, the Nile Delta, Malay Basin, Eastern Venezuelan Basin (Trinidad) and the Potwar Plateau of Pakistan 
[Law and Spencer, 1998; Burrus, 1998; Heppard et al., 1998; Powley, 1990; Nelson and Bird, 2005; Morley et al., 2011]. In 
these areas, the abnormally pressured rocks are mainly located in Tertiary and late Mesozoic sedimentary formations (146 to 
65 mya, million years ago), the depositional setting is dominantly deltaic, and the lithology is dominantly shale, as is the case 
of Block 276, where sedimentary formations, Pliocene -Pleistocene (5.3 mya to 1 1,000) is even younger. 

Additionally, Mouchet and Mitchell [1989] noted that one of major reasons of abnormal pore pressure is caused by abnormal 
formation compaction (compaction disequilibrium or under-compaction). When sediments compact normally, formation 
porosity is reduced at the same time as pore fluid is expelled. The authors added that during burial, increasing overburden 
pressure is the prime cause of fluid expulsion. They explained that if the sedimentation rate is slow, normal compaction 
occurs, i.e., equilibrium between increasing overburden and the reduction of pore fluid volume due to compaction (or ability 
to expel fluids) is maintained. They said that this normal compaction causes hydrostatic pore pressure in the formation. The 
authors explained that rapid burial, however, leads to faster expulsion of fluids in response to rapidly increasing overburden 
stress. 

Swarbrick et al. [2002] added that the overpressures generated by under-compaction in mudrock-dominated sequences may 
reveal an abnormal pore pressure change with depth. They noted that the compaction disequilibrium is known by higher than 
expected porosities at a given depth and that the porosities deviated from the normal porosity trend. 

Zhang [2011] looked at overpressure from a different angle. He plotted shale resistivity as a function of depth to explain the 
departure of pore pressure from the normal hydrostatic trend. The author identified under compaction as the main reason 
behind over pressurization. He showed that the inclined line in Fig. 4(a), above, represents the resistivity in normally 
compacted formation (normal resistivity, R n ). He noted that in the under-compacted section resistivity (R) reversal occurs, 
indicating an over pressured formation in Fig. 4(b). He showed a lower resistivity in the under compacted/over pressured 
section. In the figure, R n is the normal resistivity trendline, g v is the lithostatic or overburden stress, a e is the effective 
vertical stress, p n is the normal pore pressure, and p is the pore pressure. 


Shale resistiviiy (log scale) Shale pare pressure 



0) (b) 

Figure 4: Schematic of resistivity (a) and pore pressure (b) in an undercompacted basin, Zhang 

[ 2011 ] 
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In Fig. 5, Zhang [2011] plotted porosity as a function of depth, to prove his hypothesis. The dashed porosity profile in (a) 
represents normally compacted formation. In the over pressured section, the porosity reversal occurs (heavy line). In the over 
pressured section, porosity is larger than that in the normal compaction trendline (4> n ) and under compaction is apparent. 

In addition, Fig. 6 below shows different tracks of porosity. Porosity of shales and normal compaction trendline. Shales do 
not follow the normal compaction trend. 



Figure 5: Schematic of porosity (a) and corresponding pore pressure (b) in a sedimentary 

basin [Zhang, 2011] 


_ _ . . Pore pressure gradient (p pa) 

Porosity (u nitless) 

8 9 10 11 12 11 14 15 15 17 18 




(a) (b) 

Figure 6: Porosity and pressure prediction from the porosity [Zhang, 2011] 

In Fig. 7, Shaly PI -sand porosity profile shows the same behavior. No clear trend of compaction and porosity hovers around 
an average of 30% and stays constant with increasing depth. In Fig. 6 the normal compacted porosity should range between 
12 to 10% in PI sand. 



Figure 7: P1-sand sidewall core porosity profile 
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Permeability profiles can also be used in justifying under compaction. Fig. 8 is an illustration of permeability variation with 
depth in PI -sand. A lower permeability at the top of the formation explains the occurrence of transgressive shale, a primary 
cause of ceased compaction. Permeability is lower in the shaly PI -sand and sees a drastic shift in the bottom of the reservoir 
and indicates a facies change. 
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Figure 8: P1-sand sidewall core permeability profile 

High water saturation distribution in layer 1 (Fig. 9) is another indication of transgressive shale in large quantities and that 
PI -sand is overlain by a shale barrier, as seen in Fig. 7 and 8 above. 



FloViz 201 0.1 



Figure 9: Water saturation distribution in layer 1 


Techlog generated porosity and water saturation (Fig. 10) and volume of shale (Vshale) and net-to-gross (NTG) logs (Fig. 
11) in wells #1, #12 and W-12 also indicate that no clear trend of compaction is evident in the shaly transgressive sand and 
that the porosity is under compacted (Fig. 10). 



Figure 10: No clear trend of porosity degradation in P1-sand 
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Figure 11: Evidence of interbedded (transgressive) shale 

To better describe porosity variation with depth and under compaction, a mathematical model is used. 

IV. Mathematical Model 

The relationship between the in-situ, porosity,®, and the uncompacted porosity, ®°, is given by [Revil et al., 1998]: 

0 = 1 - (1 - 0 o )exp(^) 

where, z is a depth of interest and 


( 1 ) 


- = 9{p g -p f )<l> 0 P 


where P is the compaction coefficient with P= 4.8x1 ® psi in the Eugene Island minibasin. 

In addition, we assume that fluid pressures have not increased for reasons unrelated to compaction, (2) solid material has not 
been dissolved and removed, (3) buoyant pressures are not transmitted horizontally or vertically. 


Porosity is deduced from density logs using; 

p = (l-0)p g + 0p f (2) 

where, P is the sediment density, ^3 is the sediment grain density, and is the pore fluid density. Generally, grain density is 
taken as 165.4 lb/ft 3 and pore fluid density as 62.4. 

We assume that changes in porosity, ^0, are explicitly related to changes in effective stress, ^^eff- 

30= - 0 o P3a eff (3) 

0 R 

where, 0 is the uncompacted porosity and P is the long term compressibility. The effective stress equals the lithostatic 
stress, minus the pore pressure P: 

(T eff = P-P (4) 


Lithostatic stress is related to density of mineral grains, P 9, density of pore fluid, P f, and sediment porosity, ®> by: 

dp z = (PgC 1 " 0 ) + P /®)#Az (5) 

The interpretation method is straightforward. Equation (5) implies that instead of a decreasing porosity in line with normal 
hydrostatic compaction, a constant porosity with increasing depth implies excess fluid causing pore over pressurization. To 
make sure that over pressurization was mostly caused by under compaction; a stress analysis was conducted using COMSOL 
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to study the effect of gravity and temperature on pressure and whether the shale barrier was acting as lithostatic gradient 
absorber, limiting compaction in PI -sand. 

V. Effect of Temperature on Pressure Variation in P1-sand 


A finite element model using COMSOL multiphysics was developed to study the effect of temperature on the effective stress 
(see Fig. 12). This study will also explain whether the shale barrier overlain the sand was the main reason behind under 
compaction. A two-dimensional model of the PI -sand was built using quadrilateral elements to determine the thermal 
stresses caused by heat dissipated from the salt dome and thermal boundaries. A thermal analysis was initially conducted to 
determine the temperature distribution within the PI -sand from which thermal stresses were obtained. To mimic the 
overburden stress, a pressure gradient of 1 psi/ft was assumed and an equivalent pressure of 7500 psi was applied to the top 
of the elemental volume. Gravity of the different sedimentary strata was also considered. The eastern and western boundaries 
of the reservoir were constrained from moving laterally and allowed to move vertically; whereas, the bottom side of the 
model was constrained from moving in all directions. 


P = 7500 psi 





Figure 12: Boundary conditions used in the finite element model 


Fig. 13 shows the displacement field under mechanical and thermal loads. It can be seen that the highest vertical 
displacement of -6.43 ft occurred at the top of the section (clean sand at a depth of 7500 ft) with least displacement below the 
shale layer, which acts as a visco-plastic medium; thereby, sustaining the overburben pressure causing under compaction in 
PI -sand. Fig. 13 also shows that shale, a visco-plastic medium, absorbs most of the stress and acts as a load carrier structure 
leading to under compaction in the underlain shaly sand portion of the reservoir. 

To understand the effect of temperature on the pore pressure, the effective stress and pressure contours for both thermo - 
mechanical (gravity + temperature) and mechanical loading (gravity only) were examined. Fig.s 13(a) and (b) show the 
effective stress contour and pressure for gravity alone; while Figs. 13(c) and 13(d) portray the contours for gravity and 
temperature. A comparison between the effective stress contours reveals that temperature effect was found to contribute to 
the increase in pore pressure up to 20.5% increase throughout the whole field. 
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Figure 13: Displacement simulation contour 
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Figure 14(a): Stress distribution in simulation control volume, gravity only 
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Figure 14(b): Pressure distribution in simulation control volume, gravity only 
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Figure 14(c): Stress distribution in simulation control volume, gravity + temperature 
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Figure 14(d): Pressure distribution in simulation control volume, gravity + temperature 

A close examination of Fig.s 13(a) and (c) exhibit that Pl-sand (dark blue) is subjected to the lowest effective stress 
indicating that the pore pressure gets higher in this zone as compared to the clean sand layer (light blue); thereby, confirming 
the theory of under compaction (overpressure) presented in Equation (6). 

°ob ~ °eff + Vf ^ 

Fig. 14 shows a vertical line along which pressure values were taken at different depths. It passes through well 12 in the 
reservoir which located below the shale layer. Fig. 15 shows the pressure profile as a function of depth for both loading 
cases (mechanical only versus mechanical and thermal). The plot clearly portrays that convection is insignificant (blue line 
lays on top of the orange). Temperature generated from the salt dome has, however, affected reservoir pressure and partially 
caused over pressurization. Gray and orange lines departure of 500 in clean sand, 30 in the shale barrier, and 350 in shaly 
sand is a clear indication that over pressurization of the P-1 sand is a combined effect of under compaction and temperature. 
It is worth noting that shale layer acted as a thermal barrier, indicated by the minimal temperature effect on pressure 
variation. Moreover, the shale layer acted as a load carrier as explained by its insignificant deformation (compaction). 

L5 
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Figure 15: Schematic of well location used in COMSOL simulation 
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Figure 16: Simulated pressure profile 


VI. Conclusions 

Analysis of the Block 276 PI -sand porosity, permeability, shale content, and pressure show that departures of shale porosity 
from hydrostatic compaction trend are coherent and provide significant new information on the possible nature of under 
compaction. It was concluded that over pressuring is strongly controlled by depositional environment. A COMSOL 
Multiphysics model was also used to demonstrate that temperature effect from nearby salt diapers also contributes to over 
pressurization but not as much as under compaction, the dominant over pressurization mechanism. It was observed that the 
temperature effect on pore pressure ranges from 3 to 15%. Thereby, drilling engineers must pay attention to this increase 
whenever they plan for any drilling activity. 
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